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' Abstract. It is commonly adopted that X-rays from O stars are produced deep inside the stellar wind, and 

^SJ ' transported outwards through the bulk of the expanding matter which attenuates the radiation and affects the 

> I . shape of emission line profiles. The ability of the X-ray observatories Chandra and XMM-Newton to resolve these 

' lines spectroscopically provided a stringent test for the theory of the X-ray production. It turned out that none 

of the existing models was able to fit the observations consistently. The possible caveat of these models was 
the underlying assumption of a smooth stellar wind. Motivated by the various evidences that the stellar winds 
are in fact structured, we present a 2-D numerical model of a stochastic, inhomogeneous wind. Small parcels of 
i hot, X -ray emitting gas are permeated with cool, absorbing wind material which is compressed into thin shell 

fragments. Wind fragmentation alters the radiative transfer drastically, compared to homogeneous models of the 
same mass-loss rate. X-rays produced deep inside the wind, which would be totally absorbed in a homogeneous 
fiow, can effectively escape from a fragmented wind. The wind absorption becomes wavelength independent if 
the individual fragments are optically thick. The X-ray line profiles are flat-topped in the blue part and decline 
] steeply in the red part for the winds with short acceleration zone. For the winds where the acceleration extends 

. over significant distances, the lines can appear nearly symmetric and only slightly blueshifted, in contrast to 

' the skewed, triangular line profiles typically obtained from homogeneous wind models of high optical depth. We 

, show that profiles from a fragmented wind model can reproduce the observed line profiles from ^ Orionis. The 

■ present numerical modeling confirms the results from a previous study, where we derived analytical formulae from 
' a statistical treatment. 

I ' Key words, stars: winds, outflows - X-rays: stars -- radiative transfer 
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, 1. Introduction rates and wind velocities correctly (Pauldrach et al. 

L' ■ T . ■ 1 T However, it was pointed out early (Lucv fc Solomon [T!)70(l . 

^ • Hot massive stars possess strone; stellar wmds, as discov- j i , r • i- i. i /-/^ lu n aaa i /-> i • c 

, r / I , . . , . and later further mvestigated I (.arlberg 119801 Owocki & 

, ered by the advent of UV spectroscopy (Morton |M3. Ry^icki [», that the stationary solution for a line- 

■ One decade later it was detected that these stars emit -,. ■ j ■ ^ -u-i n j-ui- uu 

*H /Tx 1 1 „ , , rrrrwm m, • driven Wind IS unstable; small perturbation should grow 

X-rays (Harnden et al. 119791 Seward et al. 119791). Ihe on- . , r- , , rr^i- ui i i • -4. i-v. v 

- - < ) — ' ~ — * quickly and form shocks, this dc-shadowmg instability 

gin of these X-rays is debated controversially. One of the . . ,. -ii j • a -n, 

o ■'is implicitly suppressed in the LAK theory. 

first suggestions was the existence of a hot corona close to 

the stellar photosphere. However, the surrounding stellar 1-11,11 • 1 i- r ^ ■ t 

. , , T.^ 1 n 1 1 i he most detailed hydrodynamic modeling ot this line 

wind should imprint strong K-shell absorption edges on , . . , 1 1 1 1 • 1 rmvm. 

,1 ^ ^ 1.1 ^1 1 /r-i • 11- driven instability was presented by beldmeier et al. 1119971). 

the A-ray spectrum, which were not observed (Cassmelii ^, 1 1 • 1 1 • • • n n 1^ — '■ — 

n 1 n mi 1 111 j2 11 1 1 ihcse calculations show how initially small perturbations 

& b wank 119831). 1 he base-coronal model was fanaliy ruled • ,r mi • ^ • 1 

,1 r ii 1 1 r 1 i • i- 1 • • grow and form strong shocks that emit X-rays. This wind 

out because ot the lack ot characteristic coronal emission ^ ■ • n 1 1 1 1 

1. /T-, 1 T n MMn i, mi r -i 1 i shock modcl was principally able to reproduce the X- 

Imes lljaade & Lucv ll982l) . 1 heretore it became commonly ^ ^ . ^^., ,.,1 11 

1 i 1 ii i V r 1 i i 1 1 -ii • ray nux from the O supcrgiant C On A, which had been 

adopted that X-rays from hot stars are produced withm / . ,1 i- 

j_T-.li ■ 1 1^1 1 i r ii -1 • 1 • observed at that time only at low spectral resolution, 
their stellar winds, although most ot the wind gas is obvi- 
ously "cool" . 

The driving mechanism for the mass-loss from OB ^hc ultimate test for models of X-ray production in 
stars has been identified with radiation pressure on spec- ^t^ll^"^ ^^^^^ b^^^"^'^ P^^^^^^^ ^^^h the launch of the X-ray 
tral lines. A corresponding theory was developed by observatories Chandra and XMM-Newton. Their spectro- 
Castor et al. mB "CAK"). After a number of improve- g^^^P^^ ^^^^^ ^o^" ^^^^ ^'^^ "^^solve X-ray emission 
ments, this theory seems to predict the observed mass-loss 1"^^ P^o^^^^^- Surprisingly, the shape of the X-ray Imc pro- 
files turned out to be strikingly different among the few 
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6*1 OriC (Schulz et al. I^UnUjl and C OriA (Waldron 
& Cassinelli 12001(1 show profiles which appear wind- 
broadened, symmetric and unshifted with respect to 
the center wavelength; 

5 Ori A (Miller et al. I2UU2|I shows profiles which are nar- 
rower, symmetric, and also centered at the laboratory 
wavelength; 

C Pup (Cassinelli ct al. ETOTl Kahn et al. 0(M|l shows 
strongly broadened profiles which are blueshifted 
against the rest wavelength but otherwise appear to 
be symmetric. 

If X-rays are produced deep inside the stellar wind, 
they have to propagate through the absorbing cool wind 
before they can emerge towards the observer. The opac- 
ity at X-rays energies is much larger than in the visual 
and UV, due to the K-shell absorption of abundant metal 
ions and the Hen edge. Red-shifted photons of an emis- 
sion line originate in the back hemisphere of the wind, 
and suffer much stronger absorption than photons in the 
blueshifted part of the line profile. Hence, as soon as the 
absorbing wind has a significant optical depth, the line 
profiles should be blueshifted and skewed, which means 
their shape resembles that of a triangle, with the flux 
maximum at maximum blueshift (MacFarlane et al. ll991|l . 
None of the observed stars meets this prediction. 

Kramer et al. (|2003|l were able to fit eight X-ray lines 
observed in ^ Puppis with a homogeneous wind model. 
However, these fits could only be achieved by assuming 
much less absorption in the wind than expected from the 
generally adopted mass-loss rate of this star. Moreover, 
Kramer et al. H2003|l state that there is not a strong trend 
of absorption with wavelength, in contrast to what is ex- 
pected from the energy dependence of the mass absorp- 
tion coefficient. Thus even for the prototype single O star 
Puppis the wind shock model fails to reproduce observed 
X-ray emission lines consistently. At this point the ques- 
tion arose whether this model is adequate. 

We emphasize that the whole interpretation of X-ray 
spectra had been based so far on over-simplified models, 
assuming a homogeneous distribution of the absorbing ma- 
terial. By contrast, there is strong empirical evidence that 
stellar winds are in fact strongly clumped (e.g. Hamann & 
KoesterkeUnSHl Eversbcrg et al. HMSl Puis et al.^UU^. 
This supports hydrodynamic models (Owocki et al. 119881 
Feldmeier 11995(1 which show that most of the material in 
a radiatively driven wind is compressed into a series of 
dense shells by the instability. As the hydrodynamic mod- 
els are one-dimensional, they cannot tell anything about 
the lateral structure of these wind inhomogeneities. The 
de-shadowing instability acts only in the radial direction; 
on the other hand, there is no obvious mechanism which 
could synchronize the wind in different directions. The 
observed variability in wind lines is not very pronounced, 
suggesting that the dense shells break up into a large num- 
ber of fragments. 

Wind inhomogeneity alters the radiative transfer sig- 
nificantly. We have studied the effects of wind fragmenta- 



tion on the X-ray line formation in Feldmeier et al. 1(20031 
Paper I), using a statistical approach that holds in the 
limit of infinitely many fragments. We showed analyti- 
cally that the continuum wind opacity is greatly reduced 
in a structured wind, compared to a smooth wind with 
the same mass-loss rate. The line profiles we obtained 
are broad, blueshifted and flat-topped, i.e. symmetric, and 
thus of promising similarity to those observed in ^ Pup. 

In the present paper we construct a stochastic wind 
model using a numerical approach. This allows to drop 
a number of idealizations in favor of a more realistic de- 
scription. 

Statistical approach. While Paper I relied on the statisti- 
cal limit of very many fragments, we allow now for an 
arbitrary scale of fragmentation. 

Wind geometry. We study models with different spatial 
distributions of fragments, both in the radial and lat- 
eral direction. 

Wind velocity law. We abandon the restriction to a con- 
stant wind velocity and allow for so-called /3-velocity 
laws. 

Last not least, another purpose of the present paper 
is to better understand some surprising results obtained 
in Paper I. Despite of our focus on X-ray emission line 
formation, we want to emphasize their general importance 
for radiative transfer in inhomogeneous media. 

In the next section we will introduce the stochastic 
wind model. Section 3 provides the formalism for solv- 
ing the radiative transfer in a fragmented wind, as im- 
plemented in our numerical modeling. The basic effects 
of structured absorption on line profiles are illustrated 
in Sect. 4 using simplified examples. In Sect. 5 we intro- 
duce the concept of an effective opacity, and compare the 
line profiles obtained numerically from the stochastic wind 
model with the results from an analytical, statistical treat- 
ment. In Sect. 6 wc evaluate line profiles using realistic as- 
sumptions and explore the parameter space. Conclusions 
are drawn in the final Sect. 7. A forthcoming paper will be 
devoted to detailed fits of observations. 

2. The model 

2.1. Basic concepts 

Our model refers to the structure and physical condi- 
tions in the wind shock model according to the hydrody- 
namic simulations by Feldmeier et al. ((1997(1 . These dense 
fiows driven by radiation pressure on spectral lines are 
highly supersonic. The de-shadowing instability leads to 
strong gas compression. Feldmeier et al. ((1997(1 considered 
photosphcric turbulence as seed perturbation of unsta- 
ble growth. The compressed regions appear as thin shells, 
moving outward with nearly the stationary wind speed. 
Due to the lack of synchronization between different ra- 
dial directions, the shells probably break up into fragments 
of small lateral extent. The space between fragments is es- 
sentially void, but at the outer side of the dense shells exist 
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extended gas reservoirs. From their outer edge, small gas 
cloudlets are ablated and accelerated to high speed by ra- 
diation pressure. Propagating through an almost perfect 
vacuum, they catch up with the next outer shell and ram 
into it. In this collision, the gas parcels are heated and 
emit thermal X-radiation. This mechanism works within 
a limited radial range of the wind acceleration zone. 

The two structural wind components of highly com- 
pressed, cool fragments and hot X-ray emitting gas parcels 
are disjunctive. X-ray emission ceases when the wind 
reaches its terminal speed. On the other hand, despite the 
expansion due to the internal pressure, the cool fragments 
are maintained to distances of a few tens of the stellar 
radius, before they gradually dissolve into a homogeneous 
outflow. 

2.2. Assumptions 

Our stochastic model for X-ray line formation is designed 
to describe in a generic way the structures and physical 
conditions of the stellar wind, as summarized in the pre- 
ceding subsection: 

1. The flow is spherically symmetric and stationary in a 
statistical sense. The matter propagates according to 
a monotonic velocity law. 

2. The initially homogeneous wind becomes fragmented 
at some specified distance from the stellar photosphere. 
The fragmentation is maintained out to a distance 
^St^ at which the fragments dissolve into a homoge- 
neous outflow. In the radial interval {r^'^,r^^^) all 
wind material is condensed into a finite number of dis- 
crete fragments of spherical shells which propagate in 
radial direction. 

3. All shell fragments cover the same solid angle as seen 
from the center of the star, and preserve this angle 
during their propagation. In contrast to Paper I this 
solid angle is not assumed to be infinitesimally small. 

4. We consider different statistical distributions of 
the cool fragments: "concentric spheres", "broken 
spheres" , "random fragments" , and "cones" . These 
models arc described in Sect. 2.3. The statistical ap- 
proach of Paper I corresponds to the "random frag- 
ments" distribution. 

5. The X-ray emission originates from parcels of hot gas 
which arc randomly distributed. Their location is con- 
fined to the radial range r'l^,r^^. There is no self- 
absorption by emitting material and no re-emission 
of X-rays which became absorbed by cool fragments. 
Hence, emission is decoupled from absorption, which 
allows for a formal solution of the radiative transfer 
equation. 

Besides these physical assumptions, we restrict the 
model to two-dimensional symmetry. Throughout the pa- 
per the index j is used when addressing cool absorbing 
fragments, while the index i refers to parcels of X-ray 
emitting hot gas. 
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2.3. Distribution of absorbers 

We introduce cartesian coordinates as follows. The z axis 
points towards the observer, while the impact parameter 
p stands perpendicular to the latter. Alternatively, a point 
{p, z) is specified by spherical coordinates, i.e. the radius 
r — \Jp^ + z'^ and the angle — cos~^ fi between the line 
of sight and the radial vector. Hence p = z/r. 

For the velocity we adopt the standard " /3-law" 

«W = «oo (i-7)'> (1) 

where Voo is the terminal velocity, and tq is chosen such 
that v{r ~ i?*) = 0.01 Woo- R* is the stellar (photo- 
spheric) radius. For a given mass loss rate M, the continu- 
ity equation defines the density stratification via p(r) = 
(r)). We consider continuum absorption of X- 
rays by the cool wind component, caused predominantly 
by bound-free and K-shcll ionisation processes. The opac- 
ity scales linearly with the absorber density. Assuming for 
simplicity that the relative population numbers of the ab- 
sorbers do not vary with radius, the opacity x scales with 
the mass density p(r), 

= (2) 

where k is the mass absorption coefficient which can 
be calculated from the chemical composition, population 
numbers and the continuum cross sections at the fre- 
quency of the considered X-ray line. To avoid explicit ref- 
erence to these quantities, we specify the wind by its total 
radial optical depth, 

/•OC 

T* = x{r) dr . (3) 

For a homogeneous wind, the optical depth Th(p, z) 
between a point (p,z) and the observer is 

with r' ~ \/ p'^ + z'-^. 

Next we proceed to the fragmented wind. We intro- 
duce different versions of our model, according to differ- 
ent spatial arrangements of the compressed fragments (see 
Fig.CJ. 

Model "concentric spheres". First, we consider a 
model where the matter is radially compressed into Nr dis- 
crete shells forming full concentric spheres. In order that 
the time-averaged mass-fiux M of shells resembles that of 
a stationary, homogeneous wind, the probability to find a 
shell in the radius interval 7'-|-dr must scale with l/v{r). 
A corresponding set of Nr + 1 random radii is generated by 
means of the von Neumann rejection method (e.g. Press et 
al. I1992|l . All the matter located between two subsequent 
radii is compressed into a thin shell, moving with the sta- 
tionary wind speed. The mass in each shell is conserved 
with time. 
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The total mass of a homogeneous wind enclosed be- 
tween two subsequent random radii and rb is 



^-Kr^ p{r) dr — M 



dr 
v(r) 



(5) 



In the fragmented wind model, this material is swept up 
into a thin shell of infinitesimal thickness Ar located at 
Tj. The volume of this shell is Vj — An rj Ar, its matter 
density is Mj/Vj, therefore the radial optical depth across 
the shell becomes 



_rad 



Anrj 



kM 



dr 

w(r) 



(6) 



Note that the shell thickness Ar serves only as an auxiliary 
variable, and cancels in the expression for Tj. Finally, the 
location of the shell, rj G [ra,rf,], is determined by the 
condition that this radial optical depth should be the same 
as for the homogeneous wind when considering the radial 
range from where the material has been swept up: 



^rad 



kM n dr 

X(rjdr 



4tt 



v{r) 



(7) 



Requiring rJ'''' = Tj'j'"^ from Eqs. ® and 0, the shell lo- 
cation is given by 



dr 



dr 



In case of constant velocity, the shell location is just the 
harmonic mean between the boundaries, rj = y^rJTrb- 




random fragments cones 

Fig. 1. Sketch of the distribution of absorbing fragments 
for different models, for (Nr) = 1. 

Thus, if shells are distributed according to Eq. 
where r^ and rf, are two random numbers, not only the 
mass loss rate is the same as for the homogeneous wind 
but so is the radial optical depth. 

Model "broken spheres". As a next step, the absorb- 
ing spheres from the previous section arc cut into pieces 
of equal size. In the actual 2-D cross section we consider, 
we cut each ring into N§ arcs of angle 2tt/N^. Next, the 
angular positions 19 j of the arcs are randomly distributed 
between and 27r with uniform probability. Note that an 



independent set of angles -dj is generated for each previ- 
ously diced rj. The fragments partially overlap, leaving 
some radial directions free. 

Model "random fragments". In this model both coor- 
dinates rj and -dj are chosen randomly and independent 
from each other. The probability distribution is uniform 
in radius and in angle. In the statistical limit N§ — s- oo, 
and assuming constant wind speed, this model was studied 
analytically in Paper I. 

Model "cones". One might doubt whether in a real stel- 
lar wind the radial and angular distributions of fragments 
can be statistically independent, as assumed in the ran- 
dom fragment model. One expects motions to be radial 
with negligible transversal component. The inner, not yet 
fragmented zone of the wind provides a steady, homo- 
geneous, and angle-independent mass flux. The random 
fragment model conserves this mass only in the statistical 
sense, as at any moment the number of fragments in a 
given radial direction may differ from that in another di- 
rection. In a snapshot of the random fragments model, the 
radial optical depth is therefore different in different di- 
rections. In order to investigate whether this has any puz- 
zling consequences, we put forward an alternative model 
where the radial optical depth is strictly the same in all 
directions, as in the concentric spheres model. To achieve 
this, wc divide the angle interval [0, 27r] into equal 
(8) segments. Then sets of iV^ -I- 1 random radii rj are gen- 
erated independently for each of these "cones" according 
to the description given for the "concentric spheres" , and 
the fragments are placed accordingly. We emphasize that 
each fragment is confined exactly to its cone. Therefore 
the radial optical depth and the mass in each radial di- 
rection is identical to that of a homogeneous wind at any 
instant of time. 



2.4. Distribution of emission 

The emitting gas is also distributed stochastically. We as- 
sume that a large but finite number, Ncm, of small parcels 
of emitting gas are randomly distributed in the stellar 
wind in the radial interval (r""", r^^). It is assumed for 
simplicity that the emitting zones maintain their iden- 
tity during propagation with the velocity law v{r); this 
implies that the probability to find an emitting zone in 
the radial interval (r, r + dr) scales with v^^. We assume 
the same velocity law for both emission and absorption, 
and determine the location of emitting zones again by von 
Neumann's rejection method. 

After the radius r^ of an emitting parcel is selected, 
the random choice for is to be made. For a uniform 
distribution over the sphere, the probability distribution 
of must scale with sm-d. This corresponds to a uniform 
probability distribution of z in the interval (—ri,ri). As 
our model is rotationally symmetric around the z axis, 
we put the azimuth to (pi = 0. Figure|21 shows one of these 
line-emitting zones, identified by its index i, and located at 
spherical coordinates (ri,i!)i), corresponding to Cartesian 
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coordinates {zi,pi). Actually does not appear explicitly 
in our formalism, but only r.i, Zi and pi = •\/r| — zf do. 

A particular X-ray emitting zone with index i has a line 
luminosity Li. We make the following simplifying assump- 
tions. All emitting parcels of gas contain the same amount 
of matter at the same temperature. The line emission is 
powered by collisional excitation and therefore scales with 
the density-squared. During motion, each zone expands 
according to the continuity equation. Hence, the emission 
Li of zone i at radius scales with density. 



L, 



in 



rfv{ri)' 



(9) 



with a constant Lq determined by the line emissivity and 
the filling factor. 

The total X-ray line luminosity internally released by 
the wind is 



L 



X.int 



(10) 



Because of absorption by shell fragments, only a fraction 
of this radiation will emerge from the atmosphere. 



p 


















Pi 








\^>> — 




































z 



Fig. 2. Sketch of the coordinate system for the fragmented 
wind model. X-ray line photons are emitted by hot, small 
gas parcels. Shown is the emission site i located at random 
coordinates (pi, Zi), and moving with velocity v(r) in the 
radial direction. The cool wind consists of compressed and 
aligned fragments that arc also randomly distributed. The 
radiation from parcel i towards the observer at z — ^ oo has 
to cross a shell fragment j at radius Vj and covering an 



angle interval \dj 



^Ai?]. The optical depth 



across fragment j for the ray in z direction is Tj , reducing 
the line intensity 1^ by a factor exp(— Tj). 



3. Radiative transfer 

After the X-ray line emitting zones and the continuum ab- 
sorbing shell fragments have been specified, the radiative 
transfer is calculated in order to obtain the emergent X- 
ray line profile. In our model emission is decoupled from 
absorption, which allows for the formal solution of the ra- 
diative transfer equation. Formulating the radiative trans- 
fer in the usual notation (cf. Mihalas ri978|l . we define the 



emissivity r\i^ of the X-ray emitting zone with index i as 

Vi^{'^c.ml,p,Z,lf) = (11) 
S{p - p.i) 5{z - Zi) 5{p.i{ip - ifi)) (5(l/cmf - t'o) • 

In order to take advantage of the rotational symmetry 
of our model around the z axis, cylindrical coordinates 
{p,z,ip) are used here. The delta functions account for 
the negligible spatial extension of the zone and the neg- 
ligible intrinsic line width. All delta functions arc nor- 
malized with respect to the integral over their argument, 
/ 6{p — Pi)dp = 1, etc. Note that the radiation is emitted 
at the line center frequency i^q as measured in the frame 
of reference co-moving with the flow. The corresponding 
frequency ly as seen by an observer at z ^ oo is Doppler- 
shiftcd, 



J^i.obs = t'o 1 



Zi v{r,,) 



(12) 



If there were no absorption, the observer would receive an 
intensity 



I^,i{v,p, (p) 



U 
47r 



5{p~pi) 6{pi{ip - ipi)) 5(!^ - i^i,obs)(13) 



as obtained by integrating Eq. (|11|) over z. 

It may happen that X-rays emitted from the consid- 
ered parcel of hot gas have to cross one ore more absorbing 
shell fragments while traveling in the direction towards 
the observer. Thus we have first to find out which frag- 
ments cover the impact parameter pi of the emission spot. 
Imagine, the shell fragment with index j meets this con- 
dition (cf. Fig.[21l. The radial optical depth assigned to 
this fragment is Tj'^*^. If the center of this fragment has 
spherical coordinates (rj the optical depth actually 
experienced by the X-ray is 



rad 
_i 

I Mi I 



(14) 



This accounts for the longer geometrical path of rays 
which intersect the shell fragment at an inclined angle. 
For simplicity we apply a constant optical depth Tj over 
the whole shell fragment, neglecting its curvature. 

Let us compile a list J7i which contains the indices j of 
those shells which are hit by the ray from emission zone 
i. Then the total optical depth for the ray from this zone, 
Ti, is the sum over all absorption shells crossed, i.e. 



(15) 



The crucial point to understand is that the optical 
depth Ti differs essentially from the optical depth of an in- 
dividual fragment Tj . The optical depth across a fragment 
depends on the opacity x(r) and the amount of swept-up 
material. In contrast, the optical depth along a specific line 
of sight, Ti, depends in the first place on the geometrical 
distribution of absorbing fragments, which defines the list 
Ji over which the sum in Eq. (|15|1 is taken. If, for example. 
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47r2i?2 



the fragments are located in such a way that the ray from 
a particular emission zone i towards the observer does not 
cross any fragment, then = for this zone, even if all 
fragments have Tj ^ 1. Vice versa, if there are many frag- 
ments of small Tj ^ 1 on the line of sight, the summed-up 
optical depth for that ray, r^, can be significant. 

For completeness, we have to specify the absorption of 
those parts of the wind which are not considered as being 
fragmented. First, the stellar core of radius i?* is opaque, 
i.e. emission zones with Zj < and Pi < are obscured 
and produce no emergent intensity. imax 

Furthermore, we assume that the radial range between Lx = Li exp(— r^) 
R* and r"^'" is filled with homogeneous gas. Thus the ray *=i 
from an emitter at {pi , Zi ) towards the observer collects an 
additional optical depth 



this results in summing up the contributions from all emit- 
ting zones. 



Liexp(-ri) (5(i/ - i/i,obs) 



(21) 



We can calculate now the total X-ray line luminosity 
emerging from the star. For this purpose, one has to multi- 
ply the physical flux ttFi, with the surface of the reference 
sphere {AnR'^) and to integrate over all frequencies, giving 



(22) 



_int 



kM 
in 



dz 



(16) 



when traveling though this internal sphere (cf. Eq. Q). 

Here z^' = \l {rT'Y - vl and zq = max(z„ -zq). 

Moreover, for radii larger than r'^^ the fragments de- 
cayed and the wind is assumed to be homogeneous again. 
The emission from all zones which lie further inward suf- 
fers an additional absorption 



,-OUt 



with 



An 



dz 



(17) 



V ('''sh'^)'^ ~ Pi- Since we generally assume 
that is large enough for the velocity at these distances 
being constant, Eq. (|17|l becomes 



nM 
Anp 



tan - 

2 p_ 



(18) 



Therefore the optical depth along a ray between the 
emission zone i and the observer is, in general, composed 
of the sum over the crossed shell fragments, r^, plus con- 
tributions from the homogeneous parts of the wind inside 
and outside of the fragmented zone, i.e. r^"* and r™', re- 
spectively. Although this is fully implemented in our code, 
in the discussion of this paper we will refer merely to Ti 
to keep the notations simple. With t; from now to be 
understood as the total optical depth between the X-ray 
emitting zone i and the observer, the emergent intensity 
becomes 



Iuti'^,P,'P) = Ii^,i{'^,P,f) exp(-Tj) 



(19) 



The emergent astrophysical flux, at some reference ra- 
dius i?, is obtained by averaging the intensity over the disk 
of radius i?. 



1 







dp I p dif I+,^{v,p,if) 
Jo 



(20) 



We insert /j", from Eqs. and (|19|l . and perform the 
integral over the coordinates. Due to the delta functions, 



Note that the internally released X-ray line luminosity 
Lx,int from Eq. H1U|I is correctly reproduced by Eq. H21|l 
in case of zero absorption. 

Due to the discrete character of the emitting parcels 
of gas, the emergent flux is composed of delta-function 
spikes at discrete frequencies. A spectrum of finite resolu- 
tion Aiy is obtained by a corresponding binning. For each 
frequency index fc, we compile a list Tk of those emission 
zones indices i for which Vi_ohs falls into the frequency in- 
terval (i^fe — iAi/, z/fc-|-iAi/). Then, by integrating the flux 
(Eq. (|21|l ') over the frequency bin k, the binned flux at 
frequency Vk becomes 



Fii^k) 



1 



An^R'^Ai 



LiCxp(-ri) . 



(23) 



The emergent flux scales linearly with the internally re- 
leased X-ray line luminosity, Lx.int (cf. Eq. H1U|) ). It is illus- 
trative to normalize the emergent flux profllc with respect 
to the unabsorbed case. The total energy flux in the line 
for the unabsorbed case, obtained by setting all Ti ~ 
and integrating the previous equation over frequency, is 
^x.int /(47r2i?2). Normalizing the emergent flux proflle to 
this value, Eq. I|23|) becomes 

= ^ J2 L^cxp{^n) ■ (24) 



X,int 



leik 



For convenience we will display the computed line proflles 
as function of dimensionless frequency x, measured rela- 
tive to the line center and in Doppler units referring to 
the terminal wind velocity Voo-, 



(25) 



Hence the line profile is confined between normalized fre- 
quencies — 1 and 1. Using frequency bins Ax, the normal- 
ized flux flnally becomes 

^ Ljexp(-Ti) . 



F*{xk) 



Ax L 



X,int 



(26) 



If the emitting zones propagate with constant velocity and 
there is no absorption, the resulting line profile is box- 
shaped with F*{xk) = 0.5 for \xk\ < 1, so that the inte- 
gral over the whole profile correctly contains all internally 
released X-ray line luminosity (cf. Eq. I^OJ)- 
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Equation (|26|) yields the line profile for one specific 
random trial of the emitting spots and absorbing frag- 
ments. A different random trial will result in a different 
line profile, and a large number of trials (a few thousands, 
typically) is required to achieve a good statistical aver- 
age. Phenomenologically, we can consider each trial as a 
snapshot of the wind. Typical exposure times of an X- 
ray observation are much larger than the dynamical time 
scale i?*/wcx3 of the stellar wind, which corresponds to an 
average over many snapshots. Moreover, our model is 2- 
D, whereas in a 3-D wind different azimuthal angles will 
display statistically independent patterns of emission and 
absorption, whose different line profiles are also averaged 
in a flux measurement. Investigating the variability of X- 
ray lines would be interesting, but is beyond the scope of 
the present paper. 

4. Illustration of basic effects 

In this section we concentrate on simplified cases of wind 
emission and absorption to demonstrate the basic effects 
of fragmentation on line profiles. 

Homogeneous wind. The line profiles from a homoge- 
neous outflow are widely discussed in the literature (e.g. 
MacFarlane et al. ll991|l . Figurc|21 displays the line profiles 
for an optically thick, homogeneous wind of constant ve- 
locity. All emission zones are placed at the same radius, 
— Tern, which is varied for the different profiles shown in 
the figure (labels). The further out the emission is located, 
the less absorption occurs along rays in outward direction, 
i.e. in the blue part of the line profile. On the other hand, 
radiation from the back hemisphere which forms the red 
part of the line profile must always penetrate the whole 
wind. Therefore the resulting profiles appear skewed and 
blueshifted. Note that the lowest curve, corresponding to 
an emission radius of 4i?*, is zoomed 10^ times, which 
means that essentially all radiation formed deep in the 
wind is totally absorbed. 

Concentric spheres. Figure^l shows line profiles from 
the concentric spheres model, omitting the homogeneous 
inner and outer parts of the wind. Radiation emitted close 
to the stellar core is attenuated by an ensemble of con- 
centric spheres located above the emission site. The loca- 
tion of the spheres is determined by Eq. ((HJ. The optical 
depth of a shell crossed by the ray between an emitter and 
the observer. Tj, scales inversely with (Eq. H14|) ). i.e. de- 
pends on the intersection angle between ray and shell, and 
becomes large when the shell is met under a small an- 
gle. Therefore the shape of the line profile is determined 
by the range of angles under which the absorbing shell 
can be crossed by an emitted photon. For emitting spots 
placed at and absorbing fragments at r^, the smallest 
possible value of \^\ is -^/l — {I'i/rjY . Consider, for in- 
stance, the case of just one concentric shell {N,. = 1) which 
has swept up all material between r™" and r^^l^^ , moving 
with constant velocity. According to Eq. © , the shell is 
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Fig. 3. Emission line profiles from a homogeneous wind of 
constant velocity. The radial optical depth of the photo- 
sphere is = 50 for all shown profiles. All hot gas emit- 
ting the X-ray line is assumed to be located at the same 
radius rem, which is varied for the different profiles shown 
(see labels). All profiles appear skewed and blueshifted. 
The lowest curve (rem = 4i?,) is zoomed 10^ times. 
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Fig. 4. Line profiles from a wind where the absorbing ma- 
terial is compressed in TV^ concentric shells distributed 
between r^™ = 4i?, and r^'^ = 100 i?,. The emission 
originates at rem — 4i?*, the wind velocity is constant 
(v = TJoo) and the total radial optical depth across all 
shells is T, ~ 1. The homogeneous wind below 4i?* is 
omitted here. The different profiles shown correspond to 
different numbers Nr of concentric spheres (labels), whose 
randomly chosen radii are all larger than that of the emis- 
sion zone. 



located at y r™" • r""^^ in order to conserve the radial op- 
tical depth. For the example shown in Fig.Q] r™" = 4 i?» 
and r^^^ = 100 i?*, and hence the shell is located at 
Tj = 20 i?,. The emission is located at r.i = 4i?*, i.e. at a 
much smaller radius. Therefore the minimum value of \\x\ is 
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close to unity, and the optical depth is nearly the same for 
all rays crossing the shell. Consequently, in this example 
with only one shell, the line profile is nearly flat-topped, 
the normalized flux being F*{x) = 0.5 exp(— Tj) w 0.5/e. 

We increase now the number of shells iV^ into which 
the material is compressed. According to their random 
distribution, some of the shells lie closer to the minimum 
fragmentation radius than the one shell in the Nr = 1 
example. The larger Nr, the nearer is the innermost shell 
to the emitting sphere. For rays tangential to the emit- 
ting sphere, i.e. at line center, the intersection angle with 
the innermost shell becomes small. These rays experience 
strong absorption, which leads to a depression in the line 
center. 

However, these profiles are not realistic because we 
have artificially omitted any absorption inside a radius 
of 4 i?* . We now fill this region with the homogeneous in- 
ner part of the wind, while all other parameters remain as 
in Fig. 21 The effect is demonstrated in Fig.jS] If there is 
absorbing material between the stellar core and the X-ray 
emitting spherical layer, then radiation from the receding 
back hemisphere has to pass through it, whereas radiation 
from approaching part has not. Therefore, the red part of 
the line, a; > 0, is always depleted compared to the blue 
part, a; < 0. The blue part is flat-topped for a small num- 
ber of shells, but with increasing the resulting profiles 
approach the shape from a homogeneous wind. This result 
does not depend on whether the central part of the wind 
is homogeneous, or arranged in a number of concentric 
shells. 
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Fig. 5. The effect of a homogeneous inner wind. All pa- 
rameters are the same as in Fig.01 except that the cen- 
tral part of the wind between i?, and rem is now filled 
with homogeneous absorbing material. When the number 
of concentric spheres, Nr, becomes large, the emergent 
line profile converges to the line profile from a homoge- 
neous wind of the same r*. The latter is shown by the 
thick curve, which is indistinguishable from the case of 
Nr = 1000 shells. 



Optically thick broken spheres. Let us start with just 
one absorbing shell of large optical depth, assuming that 
there is no absorption in the inner part of the wind, be- 
low the emission site. When all emission is formed inside 
absorbing shell, it will be completely trapped and the ex- 
ternal observer will not register any emergent radiation. 

In the broken sphere model, we break this absorbing 
shell into N^ fragments, for simplicity all having the same 
size. In the 2-D cut each fragment is a fraction of a circle, 
i.e. an arc. These fragments are placed randomly on the 
circle, which means that the coordinate "dj of the center 
of each fragment is a random number chosen uniformly 
between and 2tt. 

Clearly, the randomly distributed fragments will over- 
lap in some directions, and leave gaps in others. Remember 
that they covered one full sphere before the latter was bro- 
ken. The radiation from inside the sphere can pass through 
these gaps entirely unattenuatcd. 

The left panel of Fig. El shows one random distribution 
of fragments. All emission is coming from just one sphere 
with radius rem, for which different values are tested (la- 
bels). Constant velocity is assumed. Therefore the cosine 
of the angle between the direction to the observer and 
to an emitting spot is the normalized frequency, fii = x, 
at which the observer receives the line photons from this 
spot. For one random distribution of the absorbing frag- 
ments, the optical depth along the line of sight from one 
specific emission spot at {pn,ri) towards the observer is 
either zero, or big. 

What is the probability that Ti{fii,ri) — 0? For solid 
spheres, the probability that a gap occurs along a line 
of sight is nil. For broken spheres, simple combinatoric 
considerations show that, on average, a fraction exp(— 1) 
of the sphere is free from absorbing fragments. 

When a photon has a probability of exp(— 1) to find a 
gap in the broken sphere, the emergent intensity is /"*" = 
/oexp(— 1), always understood in the statistical limit of 
averaging over many random settings. In other words, the 
optical depth along the line of sight is unity. Note that this 
holds for all line frequencies, because the fractional area of 
gaps in a surface does not depend on the projection angle 
under which the surface is seen. Thus the emergent line 
profile is flat-topped, independent from the radius rem of 
the emitting sphere (cf. Fig.El right panel). 

It is important to realize that the flat-topped shape is 
due to two conditions. First, the emission must be pro- 
duced in a narrow range of velocities, in order that the 
un-absorbed profile would be fiat-topped. Second, the ab- 
sorbing fragments must be thin, i.e. their radial extent 
must be small compared to their lateral size. Globular 
clouds, for instance, would not meet this condition. 

So far, we considered one single broken sphere. What if 
there are Nr concentric broken spheres, all located above 
the emission forming region? Since the i^-randomization 
of fragments is independent for each sphere, the optical 
depth for each of them is unity, and the total optical depth 
adds up to Ti = Nr- 
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Fig. 6. The broken sphere model. One opaque sphere (r* = 4) is divided into N{) fragments, which are then randomly 
redistributed in angle. One of the random trials is shown in the left panel. Constant velocity is assumed. Emission 
arises from one sphere, for which different radii rem/R* = 4, 30, 50 and 70 are tested, but the emergent line profiles 
are independent from this choice. The line profiles shown in the right panel are obtained after averaging over 10^ 
random realizations, for the case that Nj. = 1, 2 or 3 spheres were broken into fragments (labels). 



The right panel of Fig.El shows the line profiles ob- 
tained for A'^^r = 1, 2 and 3. The left panel shows one real- 
ization of a random distribution of absorbing fragments, 
whereas the line profiles are accumulated over a large num- 
ber of trials. 

The emergent flux profiles are flat-topped. Compared 
to the case without any absorption {F*{x) = 0.5), the 
flux is reduced by a factor exp(— iV^), as expected from 
our statistical considerations. The key difference between 
the line profiles in Figs.|Sl and |H1 is due to the fact that 
in the former case the transport of radiation is through 
the material, which leads to an angular-dependent opti- 
cal depth, while in the latter case the radiation emerges 
through the empty gaps, where the probability to find a 
gap does not depend on angle. 

Optically thick random fragments. So far we treated 
first the case where the radius of the absorbing (full) con- 
centric shells was chosen randomly; and second the case of 
one or more shells which were broken into fragments and 
distributed randomly in angle, while the fragments kept 
their given radial coordinate. We combine these cases now 
and distribute the individual fragments randomly and in- 
dependently in both the angular and the radial coordinate. 
For the moment the fragments are assumed to be opaque. 

The left panel of Fig.jTjshows one trial for this random 
fragments model. The total solid angle of all fragments 
corresponds to one full sphere. The absorbing fragments 
are distributed between 4.01 and 100 i?*. 

To start with, let us assume that the emitters are lo- 
cated closer to the star than any absorbing material. In 
Fig-El this corresponds to the example with rem = 4i?». 
According to the parameters of this model, there is, af- 
ter averaging over many trials, one absorbing fragment 
in each radial direction, {Nr) = 1. This is exactly the 



same situation as in the broken sphere model considered 
above. Emerging radiation suffers absorption with an opti- 
cal depth Ti ~ 1, leading to the formation of a rectangular 
line profile. 

Next let us consider emission which is permeated with 
the absorbing fragments. Again for clarity, we assume 
that the radiation is produced at a single radius rem- In 
Fig.[7| the radii of emission are chosen to be 4, 30, 50, 
and 70 i?,. When the emission takes place in outer parts 
of the wind, say at 70 i?*, there are not many fragments 
left between r^m = 70 i?* and r"^^^ = 100 i?,. As can 
be seen from the figure, there are many gaps for pho- 
tons to escape unattenuated from the front hemisphere 
towards the observer. Since we assume constant velocity 
in this example, the fragments are uniformly distributed 
in radial direction. Hence the average number of absorb- 
ing fragments along a radial ray from rem to r'^^'^ is 
(^r)(r^r'' - ^cm)/(r,T'' - ^sh"). with (Nr) = 1 in our ex- 
ample. For instance, with 70 for the radius of emission 
and 100 i?, for the outer boundary, the average number of 
absorbing fragments in radial direction is about 0.3. 

This number plays the role of an effective optical depth 
f; the probability that there is no fragment along the ra- 
dial ray from rem to the observer is exp(— f). This concept 
will be discussed on a more general basis in Sect.El See 
also the discussion after Eq. (3) in Paper I. 

For the interpretation of Fig.[71 we must understand 
how this effective optical depth depends on the value fj, 
of the emitting spot located on the sphere of radius rem 
(remember that the radiation from fi is observed at di- 
mensionless frequency x = fi). Obviously, the probability 
for a ray to hit a fragment depends only on the radial 
interval which the ray has to cross. For the front hemi- 
sphere (/X > 0) this range is r^^ — r^m, independent of 
fj,. Therefore the blue part of the profiles in Fig. is fiat- 
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Fig. 7. Random fragments model. Left panel: one random trial of the model where the absorbing shell fragments are 
randomly distributed in radius (uniformly between 4.01 and 100 i?*) and in angle. On average there is one fragment 
along each radial direction, (Nr) = 1. The fragments in this example are opaque, and the central part of the wind 
(r < 4i?*) is void. Right panel: line profiles, for emission located at rcm/R* = 4. 30, 50 and 70 (labels). For rem = 4i?*, 
all absorbers arc outside of the emission site and the resulting profile is flat-topped as in the broken sphere model (cf. 
Fig.ini). If emission is located within the absorbing wind, as in cases with rem/ R* ~ 30, 50 and 70, the red part of the 
profile is additionally absorbed in the back hemisphere, while the blue part of the line profile remains fiat-topped. 
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Fig. 8. Cones model. An absorbing sphere (i.e. a circle in the 2-D cross section) is divided into 200 equal arcs, which 
are randomly distributed in radius. In contrast to the random fragments model in Fig. 13 the fragments are kept in 
their angular position ("cone"). Otherwise, the model is the same as in Fig. 13 By eye one can not catch a big difference 
between this distribution and the one for the random fragments model. This is reflected by the line profiles shown in 
the right panel, which look very similar to those in Fig.Q The overall flux is only slightly lower than in case of random 
fragments. The biggest difference is at the blue edge of the profile. This is due to the fact that in the cones model the 
radial optical depth is preserved (here, =4), and photons emitted in almost radial direction towards the observer 
cannot find gaps to escape. 



topped, and the fiux level is cxp(— r) times the the fiux of 
the un-absorbed emission. 

If the emitters are located at rem = 4i?*, the flat- 
topped shape extends also to the red side of the profile, 
because of our somewhat artificial assumption in this ex- 



ample that the central part of the wind (r < 4i?*) is 
void. However, when the emitting spots are permeated 
with the fragments (say, at rem = 50 i?*), a ray from an 
emitting spot on the back hemisphere (/i < 0) with coor- 
dinates {pi, Zi < 0) collects optical depth from both hcmi- 
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spheres. In radius, the ray first descends from rem tor = pi 
(back hemisphere), and than chmbs from r = pi to r^^^ 
(front hemisphere) . Thus the total radial range crossed is 
+ rem — 2pi, and the optical depth is proportional 
to this range. Therefore the red part of the line becomes 
strongly absorbed and skewed. 

Cones. This model was introduced in Sect. 12. 31 Each ab- 
sorbing sphere is divided into segments (i.e. arcs in the 
two-dimensional cross section), which are randomly dis- 
tributed in radius. In contrast to the random fragments 
model, the fragments are kept at their angular position 
("cone"), and not randomized in lateral direction. This 
model is interesting, because the radial optical depth and 
the mass are conserved for each radial direction. One of 
the random trials for the distribution of absorbing frag- 
ments is shown in the left panel of Fig.|Hl As in the previ- 
ous examples, the emission is located at a unique radius 
rem, and there is no absorption below this radius. As seen 
by eye, this distribution of the fragments can hardly be 
distinguished from the one in Fig.[7| 

The resulting line profiles are shown in the right panel 
of Fig.lHl They also do not differ much from the line profiles 
obtained in the random fragment case. The overall flux is 
only slightly lower. The biggest difference is at the blue 
edge of the profiles. This is due to the fact that in the cones 
model the radial optical depth is conserved (here, =4), 
and photons emitted in almost radial direction towards 
the observer cannot find gaps to escape from absorption. 
But apart from the blue edge, the radial randomization 
of the absorbing fragments provides escape channels for 
the photons, acting in the same way as in the random 
fragments model. 

It order to verify this numerical results for 3-D, we 
made numerical experiments to illustrate the effect of ran- 
dom radial re-positioning of fragments of a sphere. We 
cover the unit sphere with spherical triangles of similar 
size, mimicking absorber fragments. Spherical coordinate 
lines are not appropriate here, due to the dependence of 
solid angle on sint9. The sphere is divided into octants. For 
each of the three grand circle segments delimiting an oc- 
tant surface, the midpoint at half arc length is connected 
to the midpoints of the other two grand circle segments. 
This divides the octant surface into four spherical trian- 
gles. For each of them this procedure (division of grand cir- 
cle segments at half arc length, and connection to nearest 
neighbors) is repeated. After six such subdivisions one ob- 
tains 4^ = 16384 triangles on the front hemisphere. These 
triangles are of similar size and shape, but not exactly 
uniform. 

Next, the radial location of these (opaque) triangles is 
chosen randomly and uniformly between r — 0.8 and 1. 
The corresponding sky projection is shown in Fig. In 
fact, the distribution appears random and uniform, except 
in the disc center and close to the limb. To quantify this, 
we divide the projected disk into 10 rings of equal radial 
extent and determine numerically the fraction of the ring 
area which is covered by at least one fragment fFie. lTTijl . 




Fig. 9. Sky projection of radially randomized fragments. 
The front hemisphere of a unit sphere is cut into 16 384 
roughly equal spherical triangles, which are subsequently 
randomly redistributed between r = 0.8 and 1.0 while 
their angular position is kept. 
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Fig. 10. Fraction of projected surface covered by radially 
randomized fragments. The sky projection of the frag- 
mented hemisphere from Fig.lHl is divided into 10 concen- 
tric rings of equal width in impact parameter p, and the 
fraction of covered (black) surface is determined numeri- 
cally. Except for the center and the limb region, the cov- 
erage is w 1 — e~^, i.e. close to the value expected if the 
fragments were randomly distributed in their angular po- 
sition, but kept on the unit sphere. 



The complement is the fractional area which remains free 
of any projected fragment, while fragments are overlap- 
ping at other locations. For central rays p <^ 1, radial 
randomization cannot move fragments out of the line of 
sight. Hence, the central region r < 0.2 remains largely ob- 
scured (small fragment overlap), r > 1. For 0.2 < r < 0.8, 
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the transmission is close to 1 — e^^, as expected if the frag- 
ments were randomly distributed laterally on the surface 
of the sphere. For r > 0.8 finally, r < 1 because frag- 
ments are radially distributed over this range. To summa- 
rize, for an observer at infinity, radial randomization in 
a spherically symmetric wind is largely indistinguishable 
from lateral randomization. 

5. Effective opacity 

An alternative approach to the stochastic numerical mod- 
eling of this paper is to introduce an effective optical depth 
fi for a random distribution of absorbing fragments by 
statistical means (cf. Paper I). Both the numerical and 
the analytical approach should yield identical results for 
the same parameters. In this section we concentrate on 
the analytical description, and compare the resulting line 
profiles with those obtained in the numerical simulations. 

Let us define the effective optical depth fj in the con- 
ventional way, as integral over an effective opacity a (di- 
mension: per length) along the path, 

adz . (27) 
The radiative transfer equation keeps its usual form 



-a ludz, 



(28) 



with the well-known solution Eq. (|19|1 . 

The effective opacity is the product of the number den- 
sity of absorbers and the (average) area of the absorbing 
fragments projected on a plane perpendicular to the line 
of sight. Moreover, we allow (in contrast to the usual defi- 
nition of atomic cross sections) for absorbers that are not 
completely opaque, but have a probability, V, that a pho- 
ton hitting the cross section is absorbed. 

The average number of shells within a unit radius 
is termed n{r). In case of constant velocity, n{r) = 
(Nr) /{rsh^~r^™) is also constant. With variable velocity, 
the radial density of shells scales with v^^ (cf. Sect. 12. 
Therefore we have 
no 



n{r) 



with 



v(r) 



no = {Nr) 




(29) 



(30) 



As can be seen from Fig.|21 the area of a fragment pro- 
jected to the plane perpendicular to the line of sight scales 
with The total area of all absorbers within a spherical 
layer of radius r and thickness dr is nir) 47rr^dr, while 
the volume itself is inr^dr. Thus the effective opacity is 

a{r,ij)^n{r) \^,\r . (31) 

Note that this opacity is not isotropic as usual, but 
angular-dependent. The probability V that a photon 
which encounters a fragment is absorbed is 
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Fig. 11. Radial optical depth of individual fragments. For 
one selected cone in the cones model, the total radial op- 
tical depth = 1 splits into Nr = 500 shells located at 
random radii rj . The optical depth rj'^^ of each fragment 
varies according to the random algorithm used to define 
the mass in each shell, and is denoted by dots. The contin- 
uous line gives the average radial optical depth fj*"^ of the 
shells, if the mass is uniformly distributed over all frag- 
ments. Further parameters of the model are r^'" = 4i?*, 
= 100 i?*, and /3 = 1 for the velocity-law exponent. 



i.e. is a function of the optical depth tj of fragments at r 
with orientation n as given by Eq. (|14() . 

As long as one avoids to integrate through the point 
where fi = 0, the path clement along the ray, dz, can be 
substituted by dr, using 



dz 



dr . 



(33) 



By inserting Eqs. (EH, lESl and ((23 into Eq. (j^Zj) we ob- 
tain for the effective optical depth between an emitting 
point at (ri,^i) and the observer 



nir^,fi,) 



n{r) (1 - e ''^) dr 



(34) 



V = 1 



(32) 



Remarkably, in contrast to Eq. ^ neither M nor k 
appears explicitly in this expression for fi. The effective 
optical depth is determined essentially by the geometrical 
distribution of fragments, as explained earlier. 

In the following we will consider the limits of optically 
thick and optically thin fragments, r, ^ 1 and tj <C 1. 

5.1. Optically thick case (tj ^ 1) 

If all cool material is assembled in a number of opti- 
cally thick fragments, the emergent radiation is only due 
to the presence of gaps between the fragment. The ra- 
diation passes through these gaps entirely unattenuated, 
but all photons which intersect a fragment are absorbed. 
Therefore the probability V is unity. In case of constant 
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velocity, the integral in Eq. 1)34(1 is taken over the constant 
n{r) and reduces to (see Eq. (4) in paper I) 



„max 



'sh 



,max 
sh 

,max 
sh 



Ti if > 

n - 2pi if fii < 0. 



(35) 



in Fig. 111! shows the run of the average optical thickness 
of fragments as function of radius (Eq. 138() ). The effec- 
tive opacity in Eq. ((37|l for the optically thin limit finally 
becomes 



The line profiles from our stochastic numerical model 
shown in Fig. 13 agree extremely well with this analytical 
result. 



5.2. Optically thin case (tj <^ 1) 

The optical depth of individual absorbers may be small 
because the total optical depth of the wind, r*, is small, 
or because the optical depth is distributed over very many 
fragments (i.e. Nr is large). If the absorbing fragments are 
optically thin, the emergent flux is both due to radiation 
which passed unattenuated through gaps, and to radiation 
which crossed fragments and was only partly absorbed. 

For small optical depths, Tj <^ 1, the exponent 
in Eq. H32|) can be expanded as cxp(— r, ) 



Consequently V w tj and Eq. H31() simplifies to 



a{r,fi) = n{r) \fi\ 



•3 ■ 



(36) 



Expressing the directional optical depth of the fragment. 



by its radial optical depth rj* 



rad 



(Eq. (HDl) yields 



a(r, /.t) = 71 (r) 



rad 



(37) 



Therefore, in the limit of optically thin fragments the ef- 
fective opacity becomes isotropic, as the usual opacity in 
homogeneous winds. 

We insert now for t'^^'^ the average optical depth of a 
fragment fj{r) at radius r. Remember that in our stochas- 
tic wind model the amount of matter swept up into an 
individual fragment is a random number. This is illus- 
trated in Fig.^3 which shows for one cone of the cones 
model how the total radial optical depth = 1 splits into 
Nr — 500 shells located at their individual radii Vj. To 
define an average optical depth, each shell shall have the 
same amount of mass. According to the geometrical dilu- 
tion, the radial optical depth of the shells scales with r~^, 
so that we can write 



12. 



(38) 



The radial integral over the effective opacity must yield 
the total optical depth t,, i.e. 



n{r) Tj(r) dr 



which leads with Eq. if^ to 
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where on the right-hand side we expressed the radial op- 
tical depth using its original definition for a homogeneous 
wind. Thus we have tq = kM / (47rno). The continuous line 



a{r, ^) 
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kM 1 
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(41) 



which is identical to the opacity x of a homogeneous wind 
given by Eq. (j2J). 

To summarize, we have demonstrated that if the op- 
tical depth across all fragments is small, tj <C 1, the ef- 
fective opacity in the fragmented wind is identical to the 
opacity of a homogeneous wind of otherwise the same pa- 
rameters. We may conclude that the homogeneous wind 
is the limiting case of a fragmented wind, when the radial 
number of fragments is big and/or the total optical depth 
is small. 
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Fig. 12. Line profiles for a model with total optical depth 
= 10, and A^^. = 10 shells in radial direction. The 
histogram-style profiles arc for the "random fragments" 
and the "cones" versions of our stochastic wind model. 
The smooth line is analytically calculated with the effec- 
tive opacity from Eq. (|31|l . and is equivalent to the statis- 
tical model in Paper I. It agrees precisely with the random 
fragment stochastic model. The more regular arrangement 
of the absorbing fragments in the cones model leads to a 
slightly lower emergent flux. In drastic contrast, a homo- 
geneous wind with the same mass-loss rate absorbs the 
X-ray line almost totally (dashed line). Further model pa- 
rameters are: constant velocity; emission distributed be- 
tween r™" — 10 i?* and r"^^ — 30 i?,; absorbing shells 
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'sh 



5.3. Comparison with the analytical statistical model 

Solving the radiative transfer equation with the effective 
opacity from Eq. H31|) is equivalent to our statistical ap- 
proach in Paper I. In this subsection we discuss how this 
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analytical treatment compares with the stochastic model 
constructed in the present paper. 

Some resulting line profiles are shown in Figs. 1121 and 
El The total optical depth is r* = 10 for both figures. 
The only difference is the number of fragments in radial 
direction, which is = 10 in Fig.^] and A^^ ~ 100 in 
Fig-El The histogram-style profiles arc for the stochastic 
wind model in the "random fragments" and "cones" ver- 
sion. The smooth line is calculated analytically using the 
effective opacity from Eq. (pTT|l . and is equivalent to that 
from the statistical approach of Paper I. It is in excellent 
agreement with the random fragment stochastic model, 
confirming its prediction of a flat-topped, blueshifted and 
nearly symmetric profile. The more regular arrangement 
of the absorbing fragments in the cones model leads to a 
slightly lower emergent flux, especially at the blue edge of 
the profile. 

For the parameters in Fig. ll2l most fragments are op- 
tically thick. Therefore the gain in emergent flux achieved 
in the stochastic wind model is dramatic. A homogeneous 
wind with the same mass-loss rate absorbs the X-ray line 
almost totally (dashed line in Fig. ll2|l . and the tiny re- 
maining line has a triangular profile. 

FigureEl shows the result for the same models, but 
with the number of fragments in radial direction being 
increased to Nr = 100. Now most of the fragments are no 
longer optically thick. Therefore the emergent line profile 
from the stochastic wind model (both versions) have a 
triangular shape like that from the homogeneous wind, 
and the gain in emergent flux is only a factor of two. 

0.004 I I I I I I I ' ' ' ' I ' ' ' ' I ' ' ' ' I I 
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Fig. 13. Same as Fig.El but with A^^ = 100 fragments in 
radial direction. Due to their larger number, most of the 
fragments are now optically thin. Therefore the profiles 
from the stochastic model are more similar to that from a 
homogeneous wind, although still stronger by a factor of 
two. 

Thus our numerical results from the stochastic wind 
model confirm fully the predictions from the statistical 
considerations in Paper I. If the wind matter is compressed 



into a small number of optically thick shells, and these 
shells break up into fragments, the effective optical depth 
of the wind is small, independent of the amount of mass 
contained in the fragments. On the contrary, a homoge- 
neous wind with the same mass-loss rate may have an arbi- 
trarily large optical depth. Therefore the gain in emergent 
intensity of an X-ray line emitted from deep inside the 
wind may be huge, comparing the fragmented with the 
homogeneous model. The transparency of the fragmented 
wind increases with decreasing number of fragments in 
radial direction. 

Wc emphasize that the line profiles from the fractured 
wind do not depend on the angular size of the shell frag- 
ments. The number of arcs in our 2-D cut of the model 
does not enter the effective opacity. The statistical limit 
of a smooth line profile can be achieved by a large N^, 
but also by a large number of random settings over which 
the results are averaged. Only for technical reasons, the 
fragments should be small enough to justify the neglect 
of curvature in our description. Moreover, in the cones 
model, too big fragments in fixed angular positions would 
spoil the frequency resolution and imprint spurious steps 
to the emergent profiles. 
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Fig. 14. Line profiles for relatively thin winds. The re- 
sult from the fragmented wind is compared to the cor- 
responding homogeneous model with the same mass-loss 
rate. The pairs arc labeled with the total optical depth 
of the homogeneous wind, t, . The stronger line of each 
pair refers to the fractured wind. The radial number of 
fragments is always (Nr) = 1. The differences between 
the fractured and the homogeneous wind become signifi- 
cant when r, ^ 1. The profiles from homogeneous models 
become weaker and skewed, while the fragmented model 
profiles remain flat-topped in the blue part. The velocity 
law exponent /3 = 1. 
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6. Parameter study 

In Sect. 21 wc employed simplified models in order to illus- 
trate and understand the efltects of wind fragmentation. 
In this section line profiles are presented for realistic wind 
parameters, calculated with our stochastic wind model in 
the "random fragments" version. 

A stellar wind is specified by a few fundamental pa- 
rameters, such as the stellar luminosity, the photospheric 
radius, the mass-loss rate, the chemical composition and 
the velocity law. Even for a homogeneous model, the cal- 
culation of the continuous opacity at a given X-ray fre- 
quency is not straightforward and requires in principle a 
Non-LTE atmospheric model. Then one can calculate the 
radial optical depth r* of the homogeneous model, as de- 
fined by Eq. (jSJ. As stated already in Sect. 12.^1 we assume 
for simplicity that the X-ray opacities scale linearly with 
density, and adopt r* of the corresponding homogeneous 
model as the free parameter of our calculations. 

The spatial structure of the inhomogeneous wind de- 
pends on the radial distribution of the wind fragments and 
of the X-ray line emitting spots. Important is the number 
of fragments in radial direction, iV^. Existing hydrody- 
namical simulations trace the fragmentation to as far as 
hundreds of stellar radii. X-ray emission is predicted to 
take place only in the wind acceleration zone, i.e. out to 
a few stellar radii. 

To keep the discussion transparent, we adopt the fol- 
lowing parameters for all figures in this section: a smooth 
wind between 1 i?* and 4 i?* , fragmentation of the wind 
between 4 i?* and 100 i?* followed by a smooth wind again 
further out; the X-ray line is emitted between 4 and 
10 i?*. The remaining parameters for the stochastic frag- 
mented wind model are the angular extension of the frag- 
ments and the number of fragments per radius. The num- 
ber of fragments in angular direction, for which we gen- 
erally use N§ — 360; is of no significance, as discussed in 
the previous section. 

Therefore, when is specified, the most influential 
parameter of the model which is not constrained by in- 
dependent means is the number of fragments per radius, 
(Nr). This parameter defines whether the fragments are 
optically thick or thin for the given mass-loss rate. 

It was shown in Section 5.2 that fragmentation be- 
comes effective in reducing opacity when fragments are 
optically thick, but cannot infiuence the radiative trans- 
fer in optically thin winds even if the number of frag- 
ments per radial direction, Nr, is small. To check the pa- 
rameter range we calculated some models with (Nr) ~ 1 
and different values of around unity. The profiles in 
Fig.El confirm that the difference between the homoge- 
neous and fragmented wind is small for the optically thin 
case T* = 0.15, and only becomes significant when r, ^ 1. 
With increasing the profiles from the homogeneous 
models become weaker and skewed, while the fragmented 
model profiles remain strong and flat-topped in the blue 
part. 
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Fig. 15. Line profiles emerging from fragmented winds 
(solid lines) and from smooth winds (dashed lines). All 
models are with r, = 1 and (N,-) = 10. The different pro- 
files refer to different velocity-law exponents (3 = 0, 0.5, 
0.8, 1, 1.5, 2, 2.5, 3, 3.5, and 4, where the higher /3 values 
yield narrower profiles. 
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Fig. 16. Line profiles emerging from an optically thick 
wind, r* = 10, compressed into one fragment per radial 
direction, (Nr) = 1. The line profiles are for different ve- 
locity laws, namely for constant velocity, and exponents 
P = 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5 (from lower to higher 
maximum) . 

The next two figures study the infiuence of different 
velocity laws on the line profile. First we consider a rel- 
atively weak wind, with optical depth = 1. The frag- 
mented wind is calculated with (Nr) = 10 shells in radial 
direction. The line profiles shown in Fig. ll5l are for dif- 
ferent values of /3 = 0.5, 0.8, 1, 1.5, 2, 2.5, 3, 3.5 and 4, 
respectively. The broadest profile is for the constant ve- 
locity case, which can be considered as the limit (3^0 
when v{r) jumps to v^o immediately at the stellar surface. 
The higher /3, the slower is the wind acceleration. CAK- 
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type hydrodynamic modeling of 0-star winds arrives at 
P ~ 0.8, typically, while certain empirical facts have been 
interpreted in favor of higher (3 values. The X-ray line pro- 
files from our stochastic wind model show some sensitivity 
to this parameter. 

At the opposite extreme, Fig.^| shows an optically 
thick wind (r* = 10) compressed into just one fragment 
per radial direction (Nr) ~ 1. The line profiles are again 
for different velocity-law exponents f3 ranging from con- 
stant velocity (/3 — > 0) to /3 = 4.5. For high values of /? 
the lines are nearly symmetric, and only slightly shifted. 
Due to the high optical depth, the emergent lines from 
the corresponding homogeneous models are far too weak 
to be plotted to the same scale. 
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Fig. 17. The gain in emergent line flux by wind fragmen- 
tation. The ratio of the integrated line flux from frag- 
mented and homogeneous wind models, TZ = F^^^/F^, is 
shown on a logarithmic scale. The curves show model se- 
ries with different optical depths (labels). Along each 
series the number of shells in radial direction, (Nr) , is var- 
ied. The gain factor TZ can be huge for very thick winds, i.e. 
high r*, but is reduced with increasing number of radial 
fragments, (Nr). All models are for a velocity-law expo- 
nent f3 = 0.8. 

The most important consequence of the wind fraction- 
ing is the drastic reduction of the wind opacity. Therefore, 
emergent X-ray line fluxes can be much higher than from 
homogeneous winds of the same mass-loss rate. The gain 
factor, i.e. the integrated emergent line flux ratio between 
a fragmented wind and the corresponding homogeneous 
wind, TZ ~ log i^f*g/i^i*, is shown in Fig.EI The fluxes are 
integrated over the whole line profile. The figure displays 
the dependence of TZ on the radial number of fragments, 
(Nr). The different series of models have different optical 
depths T, (labels). The plot demonstrates that fragmen- 
tation can increase the emergent line flux by orders of 
magnitude. The effect is extremely pronounced for very 
thick winds, i.e. high r*. Radial fragmentation into very 
many shells reduces the gain factor. 



7. Discussion and conclusions 

A main motivation for our fragmented wind model is to 
explain the observed X-ray line profiles from O stars, 
which cannot be reproduced so far. We may check now 
whether we came closer to this goal. FigurelTSl shows the 
profile of the Nex line from ( Orionis, as observed with the 
Chandra satellite (histogram). The smooth line is a ten- 
tative comparison with one of our calculations. From the 
known mass-loss rate of the star, the continuum at the fre- 
quency of this line may be optically thick (Waldron et al. 
I2001|l . Therefore we select a suitable profile (with (3 = 2) 
from the set in Fig.^lwith r, = 10. The dimensionless 
frequency is converted into wavelength according to the 
terminal wind velocity of C, Orionis (2100 kms~^), and the 
profile is convolved with a Gaussian of 0.023 A FWHM to 
account for instrumental broadening. The good agreement 
just demonstrates that our fragmented wind model can re- 
produce emergent line profiles of the observed shape, even 
if r» is large. For a systematic fit, all available line profiles 
should be considered in order to adjust the parameters of 
the model consistently. This will be subject of a forthcom- 
ing paper. 

There are two aspects in the problem of fitting the ob- 
served line profiles. One is to reproduce the line shape, 
another is to provide for the observed line flux. Ignace 
& Gayley (|2002(l considered profile shapes for optically 
thick emission lines. They showed that these lines are 
blueshifted, and slightly narrower than optically thin lines. 
The lines have a universal shape which is nearly symmet- 
ric and insensitive to the level of continuous absorption. 
Nevertheless, the flux in the line scales inversely with op- 
tical depth of the cool wind, and therefore the problem of 
missing opacity is still unresolved in their approach. The 
same problem is highlighted in Kramer et al. (|2003|l , which 
up-to-day is the most consistent attempt to fit the ob- 
served lines. Our study shows that incorporation of wind 
fragmentation can resolve the controversy. 

A burning question is, why stars with existing high- 
resolution X-ray spectra show qualitatively different pro- 
file shapes, as listed in the introduction. At present we 
suggest that only C Pup is a safe prototype for an isolated 
strong stellar wind. Indeed, speckle interferometry led to 
the discovery that 6*^ Ori C and C Ori A, are binary systems 
(Hummel et al. 20001 Weigeh et al. HM^ . and 5 Ori A is 
a well-known multiple system (see e.g. Miller et al. I2002|) . 
Moreover, 6^ OriC has a strong magnetic field in excess 
of IkG (Donati et al. I2002|l . The physical processes lead- 
ing to the emission of X-rays can differ between objects. 
Nevertheless, we claim that the effects of wind fragmen- 
tation ought to be included in the modeling of emerging 
line profiles independent of the particular emission mech- 
anism. 

The scenario chosen as a framework of the model pre- 
sented here is that of Feldmeier et al. 119971 where X-ray 
are produced in the collision of fast cloudlets with dense 
shells. The prominent feature of this hydrodynamic simu- 
lation is that the X-ray emitting plasma is always located 
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Fig. 18. Profile of the Nc x line from ( Orionis, as observed 
with the Chandra satellite (histogram) . The smooth line is 
a tentative comparison with one of our calculations, tak- 
ing the profile for = 10 and (3 = 2 from Fig. El The 
theoretical profile is scaled to the terminal wind velocity 
of C Orionis (2100 km s~^) and convolved with the instru- 
mental profile, a Gaussian of 0.023 A FWHM. The good 
fit demonstrates that our fragmented wind model can re- 
produce emergent line profiles of the observed shape. 

at the starward face of the cool absorbing fragments. It 
was already shown in Paper I, and is confirmed here, that 
this configuration leads to strong depletion of the central 
part of the line, an effect which is certainly not observed. 

The following conclusions can be drawn from our nu- 
merical modeling of the X-ray line emission from inhomo- 
geneous, fragmented stellar wind: 

1. Fragmentation drastically reduces the effective opacity 
of the wind. Therefore X-rays produced deep inside 
the wind can effectively escape, which would be totally 
absorbed in a homogeneous wind of the same mass-loss 
rate. 

2. Absorption in a fragmented wind becomes effectively 
independent of the mass absorption coefficient and 
therefore of the wavelength, if the fragments are op- 
tically thick. 

3. The line profiles from the fragmented wind model 
exhibit a variety of shapes, ranging from broad, 
blueshifted and fiat-topped to narrow and nearly sym- 
metric, and are thus of promising similarity to obser- 
vations. 

4. The possibility of a fiat-topped, i.e. almost symmetric 
blue part of the line profile is due to the model assump- 
tion that the fragments are flat, have small thickness, 
and are aligned perpendicular to the radial flow direc- 
tion. 

5. The effect of fragmentation is significant when the in- 
dividual fragments are optically thick. Therefore the 
effect depends on the average number of fragments per 
radial direction. This parameter can be empirically re- 



stricted by a detailed line fit and has interesting impli- 
cations for the theoretical understanding of the wind 
hydrodynamics . 

6. The lateral size of the wind fragments has no influence 
on the emergent lines. 

7. If the fragments are strictly confined to a pattern of ra- 
dial cones while moving outwards, which ensures strict 
mass conservation in each radial direction, the result- 
ing line profiles differ only little from the random frag- 
ment model where the absorbing shell fragments are 
randomly distributed in both radial and angular coor- 
dinates. 

8. The numerical modeling featuring stochastically ar- 
ranged emitting parcels of gas and absorbing fragments 
confirms the results of Paper I, where analytical formu- 
lae have been derived from a statistical treatment. 
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